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Abstract 

The effect of electron-nuclear spin interactions on qubit operations is investigated for a qubit 
represented by the spin of an electron localized in a self-assembled quantum dot. The localized 
electron wave function is evaluated within the atomistic tight-binding model. The magnetic field 
generated by the nuclear spins is estimated in the presence of an inhomogeneous environment 
characterized by a random nuclear spin configuration, by the dot-size distribution, by alloy disorder, 
and by interface disorder. Due to these inhomogeneities, the magnitude of the nuclear magnetic 
field varies from one qubit to another by the order of 100 G, 100 G, 10 G, and 0.1 G, respectively. 
The fluctuation of the magnetic field causes errors in exchange operations due to the inequality of 
the Zeeman splitting between two qubits. We show that the errors can be made lower than the 
quantum error threshold if an exchange energy larger than 0.1 meV is used for the operation. 
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I. INTRODUCTION 



Quantum computers (QC) hold the promise of solving problems that would otherwise be 
beyond the practical range of conventional computers. Natural candidates for the funda- 
mental building block of quantum computers (qubit) are the electronic and the nuclear spin, 
since they have a well defined Hilbert space and a relatively long decoherence time compared 
to the orbital degrees of freedom. Several QC implementations have been proposed based 
on the use of the spin degree of freedom, such as usingthe nuclear spins in a molecule 111] and 



in c^ta! latticed the nnclea, spins of donors in Sifl and in endohedra, M e re n eSli and 
the spin of electrons confined in quantum dotsN, rc| or donors. 71 While small-scale quantum 

nun 

computing has been demonstrated with a few qubits, using the nuclear spin [8, 9|, 11C| or 
trapped ions, [ill large-scale quantum computing with many qubits used in parallel is yet to 
be demonstrated. 

Solid-state spin-based QC architectures are in principle scalable to many qubits. However, 
they are intrinsically inhomogeneous due to defects, impurities, interfaces, etc. Although 
an inhomogeneous environment will cause inaccuracy in the qubit operations, fault-tolerant 
error-correction schemes will compensate for quantum errors if the error occurrence rate is 
smaller than 10~ 3 — 10~ 4 per operation. Q| For spin-based QC architectures, the single- 
qubit operation typically uses the Zeeman coupling to an external magnetic field (g/isS • B), 
while the two-qubit operation relies on the exchange interaction between two spins (JSi ■ 
S2). Therefore, if the inhomogeneous environment causes fluctuations in the local magnetic 
field, it will lead to errors in the quantum operations. For example, a recent study lj| 
has shown that a "swap" operation with a magnetic field fluctuation AB yields an error 
of (qhbAB/J) 2 . It is thus crucial to examine whether the proposed solid-state spin-based 
QC implementations are scalable within the quantum error limit, in presence of a realistic 
inhomogeneous environment. As a prototype of this examination, the present paper focuses 
on the scalability of architectures where the qubit is represented by an excess electron spin 
localized in a self-assembled InAs quantum dot. An array of self- assembled quantum dots 
is an excellent candidate for a scalable QC architecture because recent advances in the 
fabrication technology have substantially improved the control of size and location of the 
nanostructures. 

Among the many sources of inhomogeneity in the local magnetic field for an InAs quantum 
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dot embedded in a GaAs buffer, the electron-nuclear spin interaction is responsible for the 
largest magnetic field fluctuation. All the nuclei in this nanostructure possess a nonzero 
magnetic moment, and the number of nuclei interacting with the electron spin is in the 
range 10 4 — 10 6 . The effective magnetic field generated for each electron by such a large 
number of nuclear spins varies from dot to dot due to the random nuclear spin orientation, the 
dot-size distribution, alloy disorder, and interface disorder. In this paper, we first estimate 
the fluctuations in the effective nuclear magnetic field resulting from the hyperfine coupling, 
calculating the electron density using the atomistic tight-binding model and including the 
fluctuation of nuclear magnetic moments due to the inhomogeneities in the environment. The 
tight-binding model is ideally suited for the description of alloy and interface disorder with 
atomistic resolution, which enables us to study the microscopic effect of the inhomogeneous 
environment on the electron densities. In a second stage, we evaluate the effect of the 
resulting fluctuations in the nuclear magnetic field on single qubit and two-qubit operations. 

The paper is organized as follows. Section HU describes the treatment of the electron- 
nuclear spin interaction within the tight-binding model. Section IHII describes the spatial 
fluctuation of the effective nuclear magnetic field acting on the electron spin from one quan- 
tum dot to another due to the inhomogeneous environment. Section HVl discusses the effect of 
the spatial fluctuation of the nuclear magnetic field on qubit operations. Finally, Section IVl 
summarizes the results of this work. 

II. ELECTRON-NUCLEAR SPIN INTERACTION 

The electron-nuclear spin interaction originates from the coupling of a nuclear magnetic 
moment to the magnetic field generated by an electron magnetic moment (or equivalently 
from the coupling of an electron magnetic moment to the magnetic field generated by a 
nuclear magnetic moment). Both the spin and orbital angular momentum of the electron 
contribute to its magnetic moment. However, the conduction electron wave function for InAs 
self-assembled dots is mostly (>90%) composed of s-symmetry atomic orbitals and hence 
the orbital magnetic moment of the electron can be ignored. As a result, the remaining 
electron-nuclear spin interaction is described by the hyperfine Fermi contact interaction: 




(1) 



j 
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where \ib and hn are the Bohr magneton and the nuclear magneton, and gj is the g factor 
of the j'th nuclear spin. S and Ij are the spin operators for the electron and the jth nucleus, 
and r and Rj are the position vectors for the electron and the jth nucleus. Since the energy 
of the hyperfme interaction (<0.1 meV) is much smaller than the energy spacing between 
the quantized electron levels (about 10-100 meV), the hyperfme Hamiltonian for a given 
electron level can be approximated with first order perturbation theory as: 

~ l()7r ~ a 

h H f = — HBHN^gj\^(Rj)\ 2 (S -ij) 

= E^(s-i,), (2) 
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where ip(Rj) is the electron wave function at nuclear site Rj, and Aj is the effective hyperfine 
coupling constant between the electron and the jth nuclear spin. 

The coupling constant Aj is proportional to the square of the electron wave function at 
a nuclear site: 

Aj = —fiB^Ngj |^(Rj)| 2 • (3) 

Within the tight-binding model, the electron wave function is expressed as a linear combi- 
nation of atomic basis orbitals 0(r — Rj). The present tight-binding model includes sp 3 d 5 s* 
basis orbitals. |lj] Therefore, the total electron density at nuclear site Rj is given by 

|^(R,)| 2 = |^(0)+^(0)| 2 , (4) 

where aj and (3j are the tight-binding coefficients for s and s* orbitals centered at site Rj, 
respectively. In terms of the effective mass approximation, the tight-binding coefficients 
loosely speaking correspond to the envelope functions while the tight-binding orbitals cor- 
respond to the Bloch wave functions. 

The tight-binding coefficients aj and j3j depend on dot geometry, material, strain profile, 
alloy disorder, etc. The geometry of a quantum dot grown by Molecular Beam Epitaxy 
varies widely with the growth condition, [lfii llfil . \\\ Based on the experimentally achievable 
geometries, we model a lens-shaped self-assembled InAs dot with diameter 15 nm and height 
6 nm, as shown in Figure^ Since the dots are embedded in a GaAs matrix, InAs/GaAs self- 
assembled dots are strongly strained due to the large lattice mismatch of 7% between InAs 
and GaAs. The equilibrium atomic positions under strain are calculated with an atomistic 
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FIG. 1: Geometry of the self-assembled quantum dot modeled in this paper. The quantum dot is 
lens shaped with a base diameter of 15 nm and a height of 6 nm. The lines x and z are the lateral 
and vertical axes along which the spatial distributions of Aj are plotted in Figure [2j 

valence force field model Q]. The strain effect on the electronic structure is captured by 
adjusting the atomic energy levels with a linear correction that is obtained within the Lowdin 
renormalization procedure. [iJ . We also modify the nearest-neighbor coupling parameters 
for the strained structures according to the generalized version of Harrison's d 2 scaling law 
and the Slater-Koster direction-cosine rules. 20, |21[ 

In the empirical tight-binding model, (f) s (0) and S *(O) are unknown because the model 
determines the Hamiltonian matrix elements without introducing the real-space description 
of the basis orbitals. For this work, the densities of the basis orbitals at a nuclear site 
are determined empirically using measurements of the Overhauser shift of the electron spin 



resonance. 



22j The details of determining the densities are given in the Appendix. 



With the resulting aj, (3j, V's(O); and (0), Aj is calculated according to Eqs. (jHJ) and 
(HJ). The spatial distributions of Aj along the directions of the dot diameter and dot height 
are plotted in Figured The maximum value of Aj (7 neV) is found at the As nucleus located 
at the center of the quantum dot. The Aj value associated with an As nucleus is about 1.7 
times larger than that associated with the In and Ga nuclei. This large difference is due to 
the larger electron density on anions than on cations. The global distribution of Aj reflects 
the localization of the electron density. Although the electron confinement inside the dot is 
quite effective along the radial direction, along the vertical axis the electron density extends 
farther outside the dot. This causes the electron spin to interact with a large number of 
nuclei outside the dot in addition to the interaction with the nuclei inside the dot. The 
number of nuclei for which Aj is larger than 0.01 * max(A,) is about 60000, whereas the 
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FIG. 2: Spatial distributions of the hyperflne coupling coefficient (Aj) for an InAs quantum dot 
embedded in a GaAs buffer (a) along the x axis and (b) along the z axis of Fig. ^ The coupling 
coefficient Aj is given by Eq. Q and is proportional to the electron density. The dashed lines 
indicate the interface between the InAs dot and the GaAs buffer. 

number of nuclei inside the dot is about 30000. 



III. NUCLEAR MAGNETIC FIELD 

The electron-nuclear spin interaction can be expressed in terms of an effective nuclear 
magnetic field Bjy acting on the electron, which is obtained from Eq. (J2J) and is defined as 

B* = -^-£^&>- (5) 

9e^B j 

Here g e is the Lande factor of the electron and the matrix element (Ij) is taken over the 
nuclear spin state. The effective nuclear magnetic field is thus determined by the spatial 
distribution of Aj and by the nuclear spin orientation (Ij). For the quantum dot modeled 
here, we find that is of the order of 0.01 T when the nuclear spins are unpolarized, and 
of the order of 1 T when the nuclear spins are polarized. The effective nuclear magnetic field 
Bjv fluctuates spatially from dot to dot. The spatial fluctuation arises from inhomogeneities 
in the environment, such as the random nuclear spin orientation, the dot-size distribution, 
the alloy and interface disorder. 



6 



When the nuclear spins are unpolarized, the magnitude and direction of are randomly 
distributed. The fluctuation of B^r due to the random nuclear configuration is given by 



1 



9e^B 
1 



22 Aj((tj) el3s — (Ij,-)e ns ) 



£^(/, + l)u, (6) 



9e^B 

where (• • -^ns is an average over the ensemble of dots, and u is the unit vector in a random 
direction. Note that although the nuclear spin lj orientation is changing within the ensemble, 
the coupling constants Aj remain unchanged under the assumption that dot geometry and 
atomic configuration are identical for all the dots. 

The fluctuation AB^v due to the random distribution of nuclear spin orientations can 
be suppressed by polarizing the nuclear spins. However, even when the nuclear spins are 
fully polarized, the magnitude and direction of Bjv can still be broadened by other inho- 
mogeneities in the dot array. For example, an ensemble of self-assembled quantum dots has 
typically about a 10% size distribution, which is inherent to the non-equilibrium molecular 



bean epitaxy growth process. |23J] When the dot size changes, the effective number of nuclei 
interacting with the confined electron and the spatial distribution of Aj change. This leads 
to different values of B^v for dots with different sizes. The fluctuation in B^r due to the size 
distribution will be estimated here by comparing calculated values of Bjy for three different 
dot geometries with base diameter and height values of (14 nm, 5.5 nm), (15 nm, 6 nm), and 
(16 nm, 6.5 nm), respectively. From the smallest to the largest dot, the number of nuclei 
inside the dot increases from 22304 to 35161. 

We further consider two additional sources for the broadening of B^r in self-assembled 
dots: alloy and interface disorder. The alloy disorder stems from the fact that a large number 
of atomic configurations will yield the same compositional ratio in an InGaAs quantum dot. 
The interfaces between an unalloyed InAs dot and the GaAs buffer shows In-Ga intermixing 
over a length scale of 1.25 nm, which is the origin of interface disorder. |24] Alloy and interface 
disorder lead to the broadening of Btv in two ways. First, In and Ga have different nuclear 
spin quantum numbers (7i n =4.5, iGa=l-5). Second, they have different ionic potentials that 
will lead to a change in the electron densities (or Aj). 

When the nuclear spins are polarized and the dot size is uniform, the fluctuation AB^ 
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due to the alloy and interface disorder is given by 

ABa^-^-VEAW") fi> (7) 

9e^B v 

where n is the unit vector along the nuclear polarization direction, J™ is the component of 
lj along fi, and A 2 (A,/™) is the variance of A,-/™. The variance A 2 (A,/™) is studied by 
examining three different atomic configurations for an In^^Ga-rAs dot. The three atomic 
configurations are constructed by randomly choosing the cation atoms as In or Ga with 
probability 1 — x and x, respectively. The overlap between the wave functions for two 
arbitrary configurations is about 0.997, indicating that the fluctuation of the energy density 
(or Aj) is very small. Furthermore, the average of the density fluctuation per site is only 
about 10~ 4 % of the average density. Therefore, we may justifiably choose to ignore the 
fluctuation of Aj and approximate ABat as 

AB„«-^A>?A*/*n, (8) 

9e^B V 

where A 2 /™ is the variance of J™. For the case of the alloy disorder in In^^-Ga^As dots, 
A 2 /™ is calculated to be x(l — x)(Ii n — ica) 2 for all In and Ga atoms and is zero for all 
As atoms, where fi is a unit vector along the nuclear polarization direction. For the case 
of the interface disorder in InAs dots, A 2 /" is 0.25(7i n — ica) 2 for the In and Ga atoms in 
the interface region and is zero for all other atoms. Each cation atom site in the interface is 
taken to have probability 0.5 to be occupied by either an In or a Ga atom. 

After calculating the fluctuation of B^ due to the inhomogeneities in the environment 
as described above, we have obtained the following results. When the nuclear spins are 
unpolarized, a random nuclear spin configuration yields a value for ABat on the order of 
100 G. When the nuclear spins are polarized, an 10% dot-size distribution also yields AB^v 
of the order of 100 G. When the nuclear spins are polarized and the dot size is uniform, 
the alloy disorder results in fluctuations on the order of 10 G, while the interface disorder 
gives rise to fluctuations of the order of 0.1 G. These results indicate that an unpolarized 
nuclear spin configuration and quantum dot size fluctuation are the dominant sources of 
the inhomogeneous nuclear magnetic field. Since the electron localized in each quantum 
dot is immersed in a different magnetic field B^r, the Zeeman splitting (E = g e jiB^N • S) 
and precession frequency g e [j, B B N /h of each electron spin are different. This leads to a 
fluctuation in the Zeeman splitting (AE) and an ensemble dephasing with a dephasing 
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TABLE I: Nuclear magnetic field spatial fluctuation (ABn), Zeeman energy fluctuation (AE = 
g e fiB^B^), and dephasing time (T 2 * = h/g e fiBAB^) caused by various inhomogeneities for an 
InAs quantum dot embedded in a GaAs buffer. For unpolarized nuclei, each nuclear spin direction 
is chosen randomly. For dot-size fluctuations, the base diameter is set to 15±1 nm and the height 
to 6±0.5 nm. For alloy disorder, Ino.5Gao.5As dots are examined and each cation atom is randomly 
chosen to be an In or a Ga atom. For interface disorder, each cation within a 1.25 nm thick 
interface between the dot and the buffer is randomly chosen to be an In or a Ga atom, reflecting 

n 

the experimental observation of In-Ga mixing near the interface (Ref. |2J). 



Inhomogeneous Environment 



AB N (G) 



AE (eV) 



T* 2 (s) 



Unpolarized nuclei 
Dot-size fluctuation 
Alloy disorder 
Interface disorder 



100 
100 
10 
0.1 



10" b 
10~ 6 

10- 7 
10- 9 



10 -io 
10 -io 

1Q -9 
10- 7 



time T 2 * defined as h/g e fiBAB^. The values of AB^v, AE, and T 2 * resulting from the 
inhomogeneous environments studied here are summarized in Table HI 



IV. EFFECT OF NUCLEAR MAGNETIC FIELD ON QUBIT OPERATIONS 

First, we examine the effect of B^r on a two-qubit operation. The fluctuation of Btv leads 
to the inequality of the Zeeman energy in the two qubits. When a two-qubit operation such 
as a "swap" operation uses the exchange interaction JSi ■ S2, the Zeeman-energy difference 
AEz between two qubits causes an error ~ (A£ z /J) 2 .(l3| For error correction codes to be 
effective, |l2|] we require (AE/J) 2 < 10~ 4 . The Zeeman-energy fluctuation AEz due to the 
four different sources of inhomogeneity in considered in this paper ranges from 1 neV to 
1 /j,eV. Hence, the exchange energy J should be larger than 100 /xeV. At the same time, to 
prevent the electron from being excited to higher-lying orbitals, J should be smaller than the 
electron energy spacing (AE e ) between the ground and the excited orbital. The excitation 
probability due to the exchange interaction is roughly on the order of (J/AE e ) 2 . Therefore, 
J/AE e < 10~ 2 would ensure the electron to stay in the qubit space with leakage probability 
below 10- 4 . 
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The dual condition (AE Z < J<C AE e ) can be met with vertically-stacked self-assembled 
dots. 25] A recent calculation with harmonic double- well confinement potentials suggests 
that J can be varied from 10 meV to 0.1 meV as the inter-dot distance increases from 5 nm 
to 20 nm.j^ Self-assembled dots with vertical inter-dot distance as small as 2 nm can be 
easily fabricated. [3| With a given physical inter-dot distance, an effective inter-dot distance 
can be electronically tuned with gate voltages to turn on and off the exchange interaction. 
The electron energy spacing AE e of a self-assembled dot is about 50 - 100 meV, depending on 
the geometry and size of the dot. j^l Therefore, J that satisfies the dual condition is between 
0.1-1 meV, which is achievable with vertically-stacked self-assembled dots. In conclusion, 
with J between 0.1-1 meV, the error due to the inhomogeneous Zeeman energies is smaller 
than the threshold for error correction, and the qubit leakage to higher orbitals is effectively 
prevented. 

Second, we consider the effect of B N on a single-qubit operation. The single-qubit op- 
eration using the Zeeman coupling to an electron spin resonance (ESR) field (B ac cosu ac t) 
involves the tuning of the ESR field frequency to the electron-spin precession frequency or 
vice versa. This tuning can be achieved by applying a gate voltage in order to make the 
electron wave function overlap with a material having a different g factor. The tuning 
process becomes complicated in the presence of Bjv, which affects the precession frequency 
of the electron spin. The effective nuclear magnetic field B^r fluctuates in space from one 
qubit to another and evolves in time. The spatial fluctuation of Bn can be compensated 
by calibrating the gate voltage for each qubit separately. However, the temporal evolution 
of Bn is difficult to compensate since a gate calibration cannot be done immediately be- 
fore each operation. The temporal evolution is determined by many competing interactions 
such as the Zeeman coupling of the nuclear spin to the external magnetic field, the nuclear 
spin interaction with the electron spin, the nuclear spin dipolar interaction, and the nuclear 
spin-lattice interaction. Detailed studies of the temporal evolution of B^v that include these 
interactions are needed to determine how long a single-qubit gate calibration is valid. 

Here, we estimate the upper limit of the temporal change of B^v for a single-qubit gate 
calibration to be valid. We assume that a static magnetic field B n of the order of 1 T is 
applied, and that an ESR field of the order of 0.001 T is used for the spin rotation. □ □ □ 



The precession frequency of the electron spin is given by uj e = g e \iB\j (Bo + B^) 2 + (Bj^) 2 , 
where b\ and B^ are the B N component parallel and perpendicular to B , respectively. A 
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single-qubit gate will be calibrated by tuning the frequency of the ESR field uo ac to u e . After 
some time, b\ and Bjj will change by A-Bjy and Ai?^ due to nuclear spin dynamics, and 
uj e will change accordingly. This will lead to the detuning of the ESR field. The frequency 
difference u ac — u e is approximately g e /is(A.B|y + (B^/B )AB^), using that B ^> B N 
for unpolarized nuclear spins where B^ is on the order of 0.01 T. The error due to the 
detuning in the single-qubit operation is proportional to (u ac — u e ) 2 /(g e fiBB a c) 2 - To have 
the error to be smaller than the error threshold (10~ 4 ) given by the error correction code.jl^] 
Al?]y and AS^ should be smaller than 10~ 5 T and 10~ 3 T, respectively. If AB^ and AB^ 
during the time interval between the calibration and the end of QC are bigger than these 
upper limits, the single-qubit gate calibration becomes invalid, and thus QC architectures 
using only exchange interactions should be employed, j^l . I29I I30I The exchange interaction 
alone can provide universal quantum gates with a minimal cost of increasing the number of 
required physical qubits by three times and of increasing the number of gate operations by 
five to seven times. 



V. CONCLUSION 

The effect of electron-nuclear spin interactions on qubit operations is investigated here 
for a qubit represented by the electron spin localized in a self-assembled quantum dot. The 
localized electron wave function is evaluated within the atomistic tight-binding model, and 
the magnetic field generated by hyperfine coupling to the nuclear spins is estimated in the 
presence of an inhomogeneous environment characterized by random nuclear configurations, 
dot-size fluctuations, and alloy and interface disorder. Due to each of these inhomogeneities, 
the effective magnetic field felt by each electron varies on the order of 100 G, 100 G, 10 G, 
and 0.1 G, respectively. 

The inhomogeneous nuclear magnetic field causes an error in two-qubit operations due to 
the inequality of the Zeeman splitting in the two qubits. However, the errors can be made 
smaller than the quantum error threshold, as long as the exchange energy for the two-qubit 
operation is larger than 0.1 meV. Recent work indicates that an exchan ge energ y of 0.1 meV 
or larger is easily achievable with vertically stacked quantum dots.j^fil. 1.3 li Is^j At the same 
time, the large energy spacing between the ground and the excited orbital (50-100 meV) 
of the quantum dots ensures that the electron qubit stays in the ground orbital while the 
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two-qubit operation is conducted. |27[] We also present the upper limit of the temporal change 
of the nuclear magnetic field for a single-qubit gate calibration to be valid when an electron 
resonance field is used for a single-qubit operation. The changes of the nuclear magnetic 
field parallel and perpendicular to the external static magnetic field should be smaller than 
10~ 5 T and 10~ 3 T, respectively. When this condition is not met, QC architectures using only 
exchange interactions should be employed. (28, 2^| 3(| Using the exchange interaction which 



can provide universal quantum gates, quantum-dot based quantum computer architectures 
are scalable to many qubits even in the presence of inhomogeneous environments causing 
fluctuations in the effective nuclear magnetic field. 
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APPENDIX: DENSITY OF AN ATOMIC BASIS ORBITAL AT A NUCLEAR 
SITE 



In principle, a measurement of the Overhauser shift of the electron spin resonance provides 
the density of the conduction electron at the nuclear site. Although Overhauser shifts 
for bulk InAs and GaAs have not been measured, the densities of an atomic orbital at 
a nuclear site for InAs and GaAs can be deduced from the Overhauser shifts measured 
for bulk InSb.|22( The measured densities for InSb are |-?/>(0)| 2 In in lnSh ^ = 9.4 x 10 25 cm -3 
and (^(O)!^ in InSb ) = 1.6 x 10 26 cm -3 . Bulk InAs and GaAs have an ionicity similar to 
that of InSb (0.36, 0.31, and 0.32, respectively). 33J] Therefore, the distribution of the valence 
electrons between the anion and cation atoms is similar in all III-V semiconductor materials. 
This implies that the ratio of the valence electron densities in bulk is the same as the ratio 
of the defies for individual a toms: 0. FeaaH, a s a border a_a t io» we ai S o 
assume that this is true for the ratio of the conduction electron densities: 



MO) I (In in InAs) = I #j) I (In atom) 
I^WI^IninlnSb) |V>(0) atom ) ' 



(A.l) 
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in InSb) 


W°)l(Sbatom) 


Wo)l? Ga] 


n GaAs) 


IV'(O)ljGaBtam) 


l#))| 2 In] 


n InSb) 


IV<(o)| 2 Inatom) 


1^(0) |? As 


in GaAs) 


1^(0)1 (As atom) 


\m\fa, 


in InSb) 


l^(°)l(Sbatom) 



(A.2) 
(A.3) 
(A.4) 

Since the atomic density ratios are known, [35] the conduction electron densities in bulk InAs 
and GaAs can then be deduced: 

l^(0)|?inininAs)^9.4xl0 25 cm- 3 , (A.5) 

l^(0)|?AsininAs)^9.8xl0 25 cm- 3 , (A.6) 

|^(0)| 2 GainGaAs) ^5.8xl0 2 W 3 , (A.7) 

l^(0)| 2 AsinGaAs)^9.8xl0 25 cm- 3 . (A.8) 

The density |?/>(0)| 2 for each atom in bulk is related to the tight-binding orbitals S *(O) 
and <f) s (0) as follows: 

|^(O)| 2 =|«0 S (O)+/5^(O)| 2 , (A.9) 

where a and (3 are the tight-binding coefficients for the conduction-band edge wave function 
in bulk. These coefficients are listed in Table ITT1 Determining the unknown values 4> s (0) 
and 0s* (0) requires one more equation in addition to Eq. (jA.9|) . We assume that the ratio 
4> s *(0)/(j) s (0) is equal to the ratio of the corresponding atomic orbitals. The atomic orbital 
can be described by a hydrogen-like atomic orbital with an effective nuclear charge. (36, [j2| 
The ratios S *(O)/0 S (O) for In, Ga, and As atoms are 0.53, 0.44, and 0.30, respectively. 
Finally, by inserting the deduced density |-?/>(0)| 2 , the conduction-band edge coefficients a, f3, 
and the orbital ratio S *(O)/0 S (O) into Eq. (|A.9[) . the densities of the s and s* orbitals are 
obtained. The resulting densities are listed in Table ITT1 
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